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ABSTRACT 

We use a series of ray-tracing experiments to determine the magnification distribution of 
high-redshift sources by gravitational lensing. We determine empirically the relation between 
magnification and redshift, for various cosmological models. We then use this relation to estimate 
the effect of lensing on the determination of the cosmological parameters from observations of 
high-z supernovae. We found that, for supernovae at redshifts z < 1.8, the effect of lensing 
is negligible compared to the intrinsic uncertainty in the measurements. Using mock data in 
the range 1.8 < z < 8, we show that the effect of lensing can become significant. Hence, if a 
population of very-high- z supernovae was ever discovered, it would be crucial to fully understand 
the effect of lensing, before these SNe could be used to constrain cosmological models. We show 
that the distance moduli m — M for an open CDM universe and a ACDM universe are comparable 
at z > 2. Therefore if supernovae up to these redshifts were ever discovered, it is still the ones in 
the range 0.3 < z < 1 that would distinguish these two models. 



Subject headings: cosmology: theory 
— supernovae 



gravitational lensing — large-scale structure of universe 



INTRODUCTION 



High-redshift supernovae have become a major tool in modern c osmology. By measuring their appar- 
ent magnitudes, we can estimate their luminosity distances <1l (see iTonrv et al.l l2003t iBarris et al.l 12004 ; 
Riess et alJl2004l and references therein). Since the relationship between dh and the redshift z depends on 
the cosmological parameters, observati ons of distant SNe ca n constrain the cosmological model. Prior to the 
announcement of the WMAP results ( Bennett et al. 20031 ) . observations of high-z SNe provided the most 
compelling evidence of the existence of a nonzero cosmological constant. Since then, they have been used in 
combination with the WMAP data to refine the determination of the cosmological parameters. 

The luminosity distances d^ are determined by combining the observed fluxes F with estimates of 
the SNe luminosities L. Uncertainties in d^ are caused by uncertainties in L, because SNe are not perfect 
standard candles. The flux F is much easier to measure, but for distant sources the value of F might be altered 
by gravitational lensing caused by the intervening distribution of matter. For instance, a magnification fi > 1 
would result in a increase in F, and an underestimation of d^. 

Estimating the effect of lensing on the statistics of high-z supernovae is a complex problem. Using either 
an analytical model or ray-tracing simulations, we can estimate the effect of lensing of a large number of 
sources in a statistical sense. We would then need to redo the error analysis on the SNe data to include in a 
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consistent way the effect of lensing. This would be a very complex task, and in this paper we have chosen a 
much simpler approach. Our goal is not to obtain a precise estimate of the error introduced by lensing, but 
rather to assess the importance of this effect: is it dominant, important, or negligible, and for what range 
of redshift? and how does it affect the discrimination between different cosmological models? To answer 
these questions, we take at face value the published results of Type la SNe, including their error bars which 
account for every source of uncertainty but gravitational lensing. Then, we include a posteriori the effect of 
lensing, to estimate the change in the errors. This approach is not rigorous at all, and does not constitute 
a substitute for a rigorous treatment of the errors. But it has the great advantage of simplicity. We do not 
have to redo the detailed error analysis performed by the high-redshift SNe groups, and, more importantly, 
our conclusions will not be tied to any particular sample or particular data reduction and error analysis 
technique used by any particular group. We are seeking to make generic statements about the importance 
of lensing (or lack of) that are relevant to any current or future sample of high-z SNe. 

The lensing o f dista nt supernovae has been the focus of several recent studies. In an early study, 
Wambsganss et al.l (|l997l) used ray-tracing ex periments to estimate t h e effect of weak lensing on th e de- 
termination of the deceler ation parameter q . Menard fc Dalall ( 2005 ). Dodelson fc Vallinottol ( 200(1 ) . and 
Munshi fc Valagead (|2006l ) focused on SNe as a mean to study the nature of weak lensing. The issue of 
determining the cosmological p arameters for distant SNe, and how this determination is affected by lensing, 
was addr essed by Wang (2005) who used semi-analytical models to determine the magnification distribution 
function, iHolz fc Linder J 2005 ) who used Mo nte C arlo ray-tracing simu lations to study the effect of weak and 
strong lensing, and Gunnarsson et al. ( 2006h and Jonsson et al. ( 2006 ). who estimated the effect on lensing 
along individual lines of sight by considering the properties of foreground galaxies in the same direction. 
These various studie s concluded that the effect of lensing on current determinations of the cosmological 
parameters is small. lAldering et al.l (|2006l ) discussed the effect of gravitational lensing on a population of 
SNe at z > 1.7. 



What distinguishes our approach is mostly its simplicity. Our calculations depend on very few assump- 
tions, and this implies a certain amount of robustness to our results. Even though we rely on numerical 
simulations, this work should be regarded as a back-of-the-envelope calculation, whose purpose is to obtain 
a qualitative estimate of the effect of lensing on the determination of cosmological parameters by distant 
SNe. Using ray-tracing experiments, rather than a se mi-ana l ytical appr oach, enables us to ex tend our study 
to redshifts much higher than the ones considered by Wang ( 20051 ) and Holz fc Linder ( 20051 ). 



This paper is organized as follow: In §2, we describe our calculation of the magnification distribution 
P(/i), and how to estimate that distribution at any redshift z. In §3, we describe the real and mock samples 
of supernovae we use for our calculations. Results are presented in §4. In §5, we address various observational 
issues. Summary and conclusion are presented in §6. 



2. THE MAGNIFICATION DISTRIBUTION FUNCTION 

2.1. Simulations 

W e have developed a multiple lens-plane algorithm to study light propagation in inhomogeneo us uni- 
verses |PremadLMMMlAMatz^ [l9^ In 
this algorithm, the space between the observer and the sources is divided into a series of cubic boxes of 
comoving size 128 Mpc, and the matter content of each box is projected onto a plane normal to the line of 
sight. The trajectories of light rays are then computed by adding successively the deflections caused by each 
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plane. 

To use this algorithm, we need to provide a description of the matter distribution along the line of 
sight. Matter is divided into two components: background matter and galaxies. We use a P 3 M algo- 
rithm (jHocknev fe Eastwood! 1 198 if ) to simulate the distribution of background matter. The simulations 
used 64 3 equal-mass particles and a 128 3 PM grid, inside a comoving volume of size 128 Mpc. The matter 
distribution in the different cubes along the line of sight then corresponds to the state of the simulation 
at differ ent redshifts[j] We then use a Monte Carlo method for locating galax ies into the computational 
volume ( Martel. Premadi. fe Matzner 19981 Premadi. Martel. fc Matzner 1998). Galaxies are located ac- 
cording to the underlying distrib ution of backgr ound matter. Morphological types are ascribed according to 
the morphology-density relation ([Dressier! Il980l ). Galaxies are modeled as nonsingular isothermal spheres, 
with rotation velocities and core radii that vary with luminosity and morphological types. 

We consider three Cold Dark Matter (CDM) cosmological models: (1) a flat, cosmological constant 
model (ACDM) with f2 = 0.27, An = 0.73, and Hp = 71 kms'^pc" 1 . This model is in agreement with 



the results of the WMAP satellite ( Bennett et al. 20031) . (2) a low-density model with Qq = 0.3, A = 0, and 
Hq = 75kms — x Mpc _1 . (3) a matter-dominated model with f2 = 1, An = 0, and H = 75kms -1 Mpc~ 1 . For 
each model, we consider sources at 8 different redshifts: z s ~ 1, 2, 3, 4, 5, 6, 7, and 80 For each combination 
model-redshift, we performed 10-20 ray tracing experiments. Each experiment consists of propagating a 
square beam of 101 x 101 = 10,201 rays back in time from the present to redshift z s , through the matter 
distribution. The rays in the beam were widely separated, by 6 arc minutes, and therefore sampled different 
regions of space. We computed the magnification matrix A along each ray. The magnification /i is then 
given by 

M = jITT" • (1) 



det A ' 

Figure [T] shows the distribution of magnifications for the A-model. The distribution peaks at fi < 1, and is 
strongly skewed. The width of the distribution increases with the source redshift. The distributions for the 
other two models are qualitatively similar. 



2.2. Standard Deviation and Magnification Distribution 

We have determined the distributions P(/x) at some particular redshifts z s . Since SNe do not cooperate 
by going off only at these redshifts, we now want to interpolate between these distributions, to obtain P(/i) 
at any redshift. First, for each model and each source redshift z s we considered, we compute the standard 
deviation of the magnification distribution -P(/i). The values are shown in Figure [2j We use an empirical 
fit of the form 

bz 

^ = T— — , 2 

1 + cz 



where the values of b and c are given in Table 1. Thi s enables us to est i mate the values of a,, , at any 



redshift. Using the stochastic universe method (SUM) of iHolz fc Waldl (|l998l ). lHolz fc Linderl (|2005h derived 
a linear relation between cr e ff and z in the range < z < 2, for a ACDM model, where <7 e ff is the effective 



1 In practice, we combine cubes f rom different s imulations in order to avoid periodicities along the line of sight. See also the 
interesting alternative suggested by I Vale fc White! [120031) . 



The exact values of the source redshifts depend on the locations of the lens planes, which vary among models. 
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Fig. 1. — Distribution of magnifications for the ACDM model. The various curves correspond to different 
source redshifts: z s — 1 (narrowest curve), 2, 3, 4, 5, 6, 7, and 8 (widest curve). 
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Fig. 2. — Standard deviation ct m versus redshift, for all three models considered. The s olid lines show 
empir ical fits of the form — bz/(l + cz). The dotted line shows the relation derived bv iHolz fc Linder 
(|2005h . 



Table 1. Coefficients of Approximation for 



Model 


n 


A 


i? [kms _1 Mpc" 1 ] 


b 


c 


ACDM 


0.27 


0.73 


71 


0.120 


0.16 


Open 


0.30 


0.00 


75 


0.085 


0.17 


EdS 


1.00 


0.00 


75 


0.117 


0.26 
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standard deviation of a single measurement, which is not the same thing as our standard deviation a^. We 
plotted their result in Figure [5] for comparison. There is a fairly good agreement between the two methods 
at redshifts z < 2. The linear relation has a slope of 0.088. Our empirical fit for the ACDM mode{§ has a 
slope that varies from 0.120 to 0.069 in the range z = — 2. 

To determine P(fi) at any redshift z, we interpolate between the distributions we have already deter- 
mined. Consider two known distributions Pi(fJ-) and ^(a*) at redshifts z\ and z 2 that bracket z. These 
distributions satisfy the properties 

(3) 
(4) 

K)r, (5) 



Pi{n)di± = 1 , 
nPi(fj,)dfj, = 1 , 
{^-l) 2 P l {p)d i i 



where i = 1,2. We define a new distribution, 



\K) 2 - ag ]P 1 (M) + K- K) 1 ]P 2 (M) 



p( P ) = — ^ ; 7 ^ v *" 1J ^' . (6) 

T ] ' ~ '" 1 



We can easily check that this distribution also satisfies the relations ([3]H([5]). This enables us to estimate 
the magnification distribution P(p) at any redshift z. We first determine (Jfi{z) from equation and then 
substitute it in equation $Q to get at that redshift. 



THE SUPERNOVAE CATALOGS 



3.1. The Tonry et al. Sample 

Observations of high-redshift supernovae provide an estimate of the luminosity distance d^. These 
results are reported in various form in the literature. Some authors express their measurements in terms of 



effective magnitudes of distance moduli. The High-z Supernova Search Team (jTonrv et al.ll2003t I B arris et al 
2004t iRiess et al.ll2004l) express their measurements in the following form, 



log(d L H Q ) = a±S a , 



(7) 



where Hq is the Hubble constant, a is the "measurement," and S a is the "intrinsic uncertainty," which 
includes every possible source of error, except gravitational lensing. In this expression, d^Ho is expressed in 
units of kilometers per second. These authors actually use the notation (log(dffo)) for a and ± for 5 a . 



In this paper, we work with the sample of iTonrv et alj (|2003l ). This is not the most up-to-date sample, 
but it is sufficient for our purpose. This sample is comprised of 230 Type la SNe in the redshift range 
< z < 1.8, with 79 of them being located at redshifts z > 0.3 (including 5 at redshifts z > 0.9). The values 
of a and 5 a can be read directly in the 8 th and 9 th columns of their Table 8, respectively. 



3 We use n = 0.27; lHolz fc Linded ||2005|) used O = 0.28. 
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3.2. A Mock Catalog of Very-high z Supernovae 

We generated a mock catalog of 43 SNe in the range 1.8 < z < 8.1. For each "supernova," we need to 
choose a redshift z, a measured value a, and an uncertainty 5 a . There is of course no rigorous method for 
doing that, since these SNe do not exist. To provide a good coverage of the redshift range of interest, we 
use d 43 equally- s paced values of z between z = 1.8 and 8.1. To determine S a , we first plotted S a versus z, for 



the iTonrv et alJ (|2003| ) sample, to look for trends. There is a large number of SNe with z < 0.1, 6 a < 0.05. 



If we focus on the 79 SNe at redshift z > 0,3, we do not see any obvious trend, and in particular S a does not 



appear to increase with redshift. So we chose, somewhat arbitrarily, the 9 SNeQ at z > 0.828. For these SNe, 
the mean and standard deviation of the uncertainties are 5 a = 0.0631 and as = 0.0113, respectively We 
then chose the values of S a for our mock SNe randomly, by drawing them from a normal distribution with 
mean S a and standard deviation erg. This ensures a smooth transition between the real and mock samples. 

To determine a, we assume that the underlying cosmology corresponds to a ACDM universe (as sup- 
ported by the real sample). We then use 

a = cia + Aa , (8) 

where cia = log(dA-ffo) is t ne actual value of \og(dr 1 H ) in a ACDM universe, and Aa is a random number 
drawn from a normal distribution with mean and standard deviation given by the value of 5 a we just 
calculated. 



4. THE EFFECT OF LENSING ON STATISTICS OF HIGH-Z SUPERNOVAE 

4.1. Compounding the Errors 

As we explained in §1, our goal is not to perform a rigorous error analysis of the uncertainties resulting 
from the possibility of lensing, but rather to estimate a posteriori the effect of lensing on the uncertainties 
already present in the analysis. 

We estimate the effect of lensing as follows: the distances of high-z supernovae are reported in the 
literature in the format given by equation ([7]), where S a is the intrinsic uncertainty (i.e. not caused by 
lensing). The distance is related to the luminosity L and flux F by 

We use equations j7]) and ([9]) to eliminate cIl, and get 

L 1/2 H = 10 a 10 ±5 »(47rF) 1/2 . (10) 

The effect of lensing will be to modify the flux F. To account for it, we replace F by F±AF in equation (I10p . 
and perform a Taylor expansion to first order in AF, 

L^ 2 H = 10°10 ±tf -(47rF) 1/a ^1 ± ■ (11) 



4 We deliberately avoided sn97G and sn76cl, whose uncertainties are much larger than those of any SNe at comparable 
redshift. 
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This expression reduces to 



log(diiJo) = a±5 a ± 



AF 



2FlnlO 



(12) 



The last term represents the uncertainty due to lensing. For a given supernova with magnification /i, 
AF/F = fi — 1. Of course, we will never know the value of /j, for a single source. But for a large number of 
sources, we can use statistics. First, the simplest, lowest-order approximation for a "typical" value of ji is 
fi = (/i) ± (Tju = 1 ± <t m , or equivalently AF/F « er M . Equation (112p reduces to 



log{d L H a ) = a±S a ±5 fl , 



(13) 



where 6 fl (z) — a /M (z)/2 In 10 can be computed using the empirical relations plotted in Figure [2j We use the 
values of a and 5 a reported bv iTonrv et alJ (|2003l ) (their Table 8). In Figure 02 we plot the ratio S^/Sa 
versus z (left of the dotted lines). This quantity increases with redshift, but never gets higher than 0.5 for 
the Tonry et al. sample. Furthermore, we shall assume that 8 a and 5^ are statistically independent, and 
combine them in quadrature, using 

S=(S 2 a + 6ty/\ (14) 



where 5 is the total error. The contribution of lensing to this error is then of order 25% at most. 

For each supernova, we compute the quantity A (to — M) (deviation of the difference between apparent 
and absolute magnitude, relative to an empty universe), using 



A(m - M) = 5\og(d L H ) - 5log(d L H Q ) ompty = 5log{d L H ) - 5 log 



cz[l + - 



(15) 



We then average the quantities A(m — M) and S in redshift bins, using 



[Aim-M)}^ = T, iWi A{m-M)l^ iWi , (16) 
Sj = (l/^w t ) 1/2 , (17) 



where 



wt = 1/Sf , (18) 

and the sums are over all data points i in bin j (note: eq. [T7] comes from 1/6? = X^JV^f])- Notice that this 
method of averaging is much f ancier than what appears to be done in the supernovae papers. For instance, 
Figure 9 of ITonrv et all |2003[ ) shows an averaging over redshift bins which is based on the median of the 



data and apparently does not take into account the uncertainties S a on the individual supernovae. 

Figure |4] shows a Hubble diagram [deviation A(to — M) versus redshift]. The data points and error 
bars on the left hand side of the dotted lines correspond to the values given by equations (|16p and (|17[) . 
respectively. The three curves, from top to bottom, show the exact results for the ACDM, low-density, 
and matter-dominated models, respectively. The results support the flat ACDM model and exclude the 
other models considered. The other panels of Figured] show the effect of lensing (the three models have to 
be plotted separately, because the correction due to lensing, which uses the relations plotted in Fig. [21 is 
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Fig. 3.— Ratio S^/Sa versus redshift. The dotted lines separate the real data of Tonry et al. (left side) from 
the mock, high-redshift data (right side). 
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Fig. 4. — Hubble diagram showing the magnitude deviatio n Afro — M) relativ e to an empty universe, for 
the three models considered. The dotted lines separate the Tonrv et al. ( 2003h data (on the left) from the 
mock data (on the right). In the top panel, the three curves, from top to bottom, show the analytical result 
for the cosmological models (f2o,Ao)=(0.27,0.73), (0.3,0.0), and (1.0,0.0), respectively. The error bars show 
90% confidence level. The last three panels reproduce the data of the top panel, but have been corrected 
to account for lensing. Since this correction is model-dependent, the three models are plotted on separate 
panels. 
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model-dependent) . This effect is totaiiy negligible. The largest correction to the error bars is about 10% for 
the highest redshift bin, for the ACDM model. 

Clearly, the potential error introduced by lensi ng is negligible in compa rison to the intrinsic error in the 
measurement, at least for SNe at redshifts z < 1.8. iGunnarsson et al.l (|2006h and lJonsson et al.l ([20061 ) reach 
the same conclusion, We now estimate the effect of lensing on a yet-undiscovered population of very-high- z 
SNe, using our mock catalog. The ratios <5 M /<5 a are plotted in Figure O on the right hand side of the dotted 
lines. The effect of lensing rapidly becomes important at redshift z > 2, especially for the ACDM models. We 
find many SNe with |<5 M | > 8 a , that is, the correction due to lensing is larger than the intrinsic uncertainty. 

The points located on the right hand side of the dotted lines in Figure [4] shows the results for the mock 
data. The error bars get significantly bigger when lensing is included. Furthermore, at redshift z ~ 3, it 
becomes very difficult to distinguish the open, low density model from the cosmological constant model, 
because the theoretical curves intersect. Keeping in mind the caveat that the mock catalog was built under 
the assumption that the underlying cosmology was ACDM, we see that the Einstcin-de Sitter model is totally 
ruled out by SNe at z > 2, but the open model is not. Indeed, it is clear that SNe at z > 2 would be rather 
useless in distinguishing an open CDM and a ACDM model: the theoretical curves get closer, whereas the 
error bars become larger. It is, interestingly, the SNe in the redshift range 0.3 < z < 1 that would still provide 
the best discriminant between these two models, and data in that redshift range are already available. 



4.2. Monte Carlo Approach 

The calculation presented in the previous section relies entirely in the standard deviation for esti- 
mating the uncertainties caused by lensing. This approach would probably be sufficient if the distributions 
of magnifications P(/jl) were gaussian. However, for sources at large redshifts, P(fJ-) is strongly skewed, as 
Figure 1 shows. The large majority of sources are demagnified, as the light reaching the observer travels 
mostly through underdense regions of the universe, while a few sources are strongly magnified, especially 
those which happen to be aligned with a massive galaxy at intermediate redshift. 

To account for the distribution of magnifications, we consider all SNe at redshifts z > 0.9 (5 from the 
Tonry sample, 43 from the mock catalog). For each one, we determine the distribution P(fi) at its redshift 
z, using equation ([6|), and then choose a magnification jj, by drawing it randomly from the distribution P(ju). 
We then compute 5^ = AF/2F lnlO = (/i — l)/21nl0. The resulting ratios |5 M |/5 a are plotted in Figure[5] 
Comparing with Figure [3j we find only a few SNe for which this ratio exceeds unity. For all redshifts and 
all models, we find that the distributions peak at a value ^ poa k < 1 such that |^t pea k — 1| < Hence, 
setting AF/F = (instead of AF/F = fi — 1), as we did in §4.1, overestimates the effect of lensing for 
most SNe. However, the distributions are very skewed, and as a result a few SNe are highly magnified, as 
Figure [5] shows. 

Figure [5] shows the resulting Hubble diagram. Comparing with Figure [¥J we find that the error bars 
are significantly smaller. The effect of lensing is less important when we use the actual distribution of 
magnifications P(/J-), and not only its standard deviation. However, the results at z > 1.8, which assume an 
underlying ACDM cosmology, still cannot rule out the open CDM model; the error bars are still too large. 
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Fig. 5. — Same as Figure [3l but with 8^ computed using the Monte Carlo approach for all SNe (real or 
mocked) at redshift z > 0.9. 
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5. DETECTION LIMIT AND BIASING 

So far, wc have assumed that any supernova, with any value of z and fi, can be observed. This assumption 
is probably valid over the range < z < 2, which includes all the current observations. But as the redshift 
gets higher, it becomes increasingly difficult to observe supernovae with current and even future telescopes, 
because of the combined effect of the flux reduction and the light being shifted to the near infrared, ft 
might just be impossib le to detect a supernova at redshift z > 2 unless, of course, it is magnified by lensing. 
Marri fe Ferrara ( 19981 ) argue that, with magnification taken into account, it might be possible to observe 



Type II SNe at redshifts up to z — 10. Of course, if only the highly-magnified SNe are observable, this 
introduces biasing, an effect that we must take into account. 



Here, we do not want to perform a detailed analysis similar to the one of iMarri fe Ferraral (I1998T). but 



j ust to get a rough estimate of the importance of biasing. The key results of the analysis of 



Marri fe Ferrara 



(|1998l ) are shown in their Figure 7, which shows, as a function of redshift, the AB apparent magnitude in 
the J, K, L, and M bands, with and without lensing. They also plot the expected flux limit of JWSlH. 
From this figure, we see that without lensing, the apparent magnitudes are below the detection limit at high 
redshift. When lensing is included, the magnitudes in the J and K bands are above the detection limit. Their 
analysis was for Type II SNe. Type la SNe are typically 1.5 magnitudes bright er, and therefore much e asier 



to detect. Neglecting the details of the spectra, we can simply take Figure 7 of IMarri fc Ferraral (|1998| ) and 
shift all the curves upward by 1.5 magnitudes. We find that, without lensing, Type la SNe would be visible 
in the J and K band, and almost visible in the L band, at redshifts up to z ~ 8. With lensing, most SNe 
would be demagnified, but the reduction in flux is typically of the order 10 — 20%, that is, a correction of 
0.103 — 0.198 magnitudes. Hence, all high-z Type la SNe should be detectable, using the proper telescope, 
and therefore we were justified to ignore any biasing effect. This being said, the identification of Type la 
SNe requires that we obtain a spectrum, and this could be quite challenging at these extreme redshifts. 

We have assumed that the typical intrinsic uncertainties S a do not grow with redshift for z > 1.8, based 
on the absence of obvious trend at z < 1.8. If the uncertainties do grow with redshift, our conclusion that 
understanding the effect of lensing at high redshift would be weakened, in the sense that these data would 
be rather useless no matter how well lensing is understood. Furthermore, it would reinforce our conclusion 
that SNe at intermediate redshifts are more useful to discriminate between different cosmological models. 

We should also ask whether Type la SNe at redshift z = 8 can actually exist. For a ACDM model 
with a Hubble constant of 71 kms _1 Mpc _1 , this redshift corresponds to an age of the universe of 650 Myr. 
Subtracting the formation and evolutionary time of the progenitor, we are getting embarrassingly close to 
the big bang. Type II SNe would be a far better candidate for high-z SNe, since the evolutionary time of 
their progenitors are much shorter. But then only the few that are magnified would be detectable, and their 
number might be too small to perform any meaningful statistics on them. 



6. SUMMARY AND CONCLUSION 

We have performed a series of ray-tracing experiments using a multiple lens-plane algorithm. We 
have determined the distributions of magnifications P([i) for sources in the redshift range < z < 8, for 
three different cosmological models. We have used these distributions to estimate the effect of gravitational 



5 They called it NGST back then. 
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lensing on the determination of the cosmological parameters with high-redshift Type la supernovae. We 
used a generic, a posteriori approach which is not tied to any particular sample. 

We found that errors introduced by lensing are unimportant for SNe with redshift z < 1.8. These errors 
are negligible compared to the intrinsic errors already present in the supernovae data. Since those intrinsic 
errors do not prevent us from determining the cosmological parame ters, the additional er rors introduced by 
lensing have no consequences. A similar conclusion was reached by Aldering et al. ( 2006h . 



Using a mock catalog of high-z SNe, extending to z — 8.1, we showed that the effect of lensing on a 
hypothetical population of SNe at redshifts z > 2 could be quite significant, and must be understood before 
such SNe could be used to constrain cosmological models. Furthermore, the open CDM and ACDM are 
difficult to distinguish at that redshift. We showed that, even if SNe at redshift z ~ 8 were ever discovered, 
it is the SNe in the range z = 0.3 — 1 that would still provide the best discriminant between these two models. 
The data at that redshift already exist, and they support the ACDM model. 



This work benefited from stimulating discussions with Gilbert Holder, Daniel Holz, Eric Linder, Massimo 
Meneghetti, and Christopher Vale. The calculations were performed at the Texas Advanced Computing 
Center, University of Texas, and the Laboratoire d'Astrophysique Numerique, Universite Laval. HM thanks 
the Canada Research Chair program and NSERC for support. PP thanks Uros Seljak for various advice and 
hospitality during the fruitful visit at the Abdus Salam ICTP in Trieste. 



- 15 - 



REFERENCES 

Aldering, G., Kim, A. G., Kowalski, M., Linder, E. V., & Perlmutter, S. 2006, Astropart. Phys., 27, 213 

Barris, B. J. et al. 2004, ApJ, 602, 571 

Bennett, C. L. et al. 2003, ApJS, 148, 1 

Dodelson, S., & Vallinotto, A. 2006, PRD, 74, 063515 

Dressier, A. 1980, ApJ, 236, 351 

Gunnarsson, C, Dahlcn, T., Goobar, C, Jonsson, J, & Mortsell, E. 2006, ApJ, 640, 417 

Hockney, R. W., & Eastwood, J. W. 1981, Computer Simulation using Particles (New York: McGraw Hill). 

Holz, D. E., & Linder, E. V. 2005, ApJ, 631, 678 

Holz, D. E., & Wald, R. M. 1998, Phy. Rev. D., 58, 063501 

Jonsson, J., Dahlcn, T., Goobar, A., Gunnarsson, C, Mortsell, E., & Lee, K. 2006, ApJ, 639, 991 

Marri, S., & Ferrara, A. 1998; ApJ, 509, 43 

Martel, H., Premadi, P., & Matzner, R. 1998, ApJ, 497, 512 

Martel, H., Premadi, P., & Matzner, R. 2000, ApJ, 537, 28 

Menard, B., & Dalai, N. 2005, MNRAS, 358, 101 



Munshi, D., & Valageas, P. 2006, preprint (astro-ph/0601683) 

Premadi, P., Martel, H., & Matzner, R. 1998, ApJ, 493, 10 

Premadi, P., Martel, H., Matzner, R., & Futamase, T. 2001a, ApJS, 135, 7 

Premadi, P., Martel, H., Matzner, R., & Futamase, T. 2001b, Pub.A.S.Aus., 18, 201 

Riess, A. G. et al. 2004, ApJ, 607, 665 

Tonry, J. L. et al. 2003, ApJ, 594, 1 

Vale, C, & White, M. 2003, ApJ, 592, 699 

Wambsganss, J., Cen, R., Xu, G., & Ostriker, J. P. 1997, ApJ, 475, L81 
Wang, Y. 2005, JCAP, 503, 5 



This preprint was prepared with the AAS IATjrjX macros v5.2. 



